Micron-scale droplet deposition on a hydrophobic surface using a retreating syringe 
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Droplet deposition onto a hydrophobic surface is studied experimentally and numerically. A wide 
range of droplet sizes can result from the same syringe, depending strongly on the needle retraction 
speed. Three regimes are identified according to the motion of the contact line. In Region I, at 
slow retraction speeds, the contact line expands and large droplets can be achieved. In Region II, 
at moderate needle speeds, a quasi-cylindrical liquid bridge forms resulting in drops approximately 
the size of the needle. Finally, at high speeds (Region III), the contact line retracts and droplets 
much smaller than the syringe diameter are observed. Scaling arguments are presented identifying 
the dominant mechanisms in each regime. Results from nonlinear numerical simulations agree well 
with the experiments, although the accuracy of the predictions is limited by inadequate models for 
the behavior of the dynamic contact angle. 

PACS numbers: 47.55.D-,47.55.np,47.55.nk 



Contact dispensing methods of fluids are widely used 
in a variety of applications including direct scanning 
probe lithography 1] , micromachined fountain-pen tech- 
niques [1, y| , electrowetting-assisted drop deposition Q 
and biofluid dispensing applications [B|, Q . The process 
is, at first glance, straightforward and is initiated by the 
formation of a liquid bridge between the substrate and 
a dispensing syringe. As the syringe retreats, the liquid 
bridge stretches, grows and breaks, leaving a drop on the 
substrate. A seemingly simple question can be asked - 
how does the drop size depend on the syringe geometry, 
speed and the fluid properties? A comprehensive answer 
must consider the stability of the liquid bridge and the 
physics of the moving contact line at the liquid-air-solid 
interface - both difficult problems. Theoretical studies of 
liquid bridge stability date back to Rayleigh Q , and have 
been extended to include gravity and non-cylindrical ge- 
ometries [alSj- I n addition, the nonlinear dynamics have 
been solved numerically, using both 2-D (axisymmetric) 
and 1-D (slender-jet) [H |12| models. Previous work 
has concentrated on geometries in which the contact line 
is pinned at both ends of the liquid bridge [12, EH, 
and there are only a few results that couple the liquid 
bridge with a moving contact line [3, [HI- A possible 
reason for this is the difficulty in solving the flow near 
the contact line where the continuum equations are in- 
valid [lit \vh and a microscopic description must be im- 
posed (e.g. [l8(). In this letter, we focus on the physics 
of drop dispensing on a flat, smooth, hydrophobic sub- 
strate in which the contact line is free to move and is 
inherently coupled with the liquid bridge stability. Ex- 
periments and numerical simulations are used to identify 
a range of complex flow phenomena which enable the de- 
posited drop size to vary by two orders of magnitude as 
the syringe retraction speed is changed. 

In our experiment, a stainless steel syringe (typical ra- 
dius, R = 200/xm) is mounted vertically on a computer- 
controlled stage. The syringe is connected by a small 
tube to a lOcc barrel mounted on the same stage. This 



configuration maintains a constant hydrostatic head, H, 
at the syringe tip (H ~ 4cm). The fluid (a 85-15 
mixture by volume of glycerol and water) has viscosity 
fi — 84 cP and surface tension 7 = 0.063 N/m. The fluid 
exhibits a static contact angle of ~ 90° with the sub- 
strate, a smooth glass slide coated with a monolayer of 
octadecyltrichlorosilane (OTS). The syringe is brought 
down towards the substrate, stopping ~ 40/zm above the 
surface so that the meniscus touches the substrate and 
spreads, partially wetting the surface to form a stable 
drum-shaped liquid bridge (Fig. [T]). As the syringe re- 
tracts at a constant speed, U, the liquid bridge elongates 
and evolves due to the changing height, h, fluid flowing 
into the bridge through the syringe, (characterized by an 
inflow velocity, Uf, and the motion of the contact line 
between the bridge and the substrate (characterized by 
a contact line position, r, and speed, u c ). At a criti- 
cal height, h pi the liquid bridge becomes unstable and 
pinches off rapidly, leaving a drop on the surface. A high 
speed camera (Photron APX) equipped with a 5X Mi- 
tutuyo lens was used to capture the motion of drop dis- 
pensing at frames rates up to lOkfps, with a resolution 
of 3.33//77i/pixel. The experiment was carried out using 
several syringe diameters, hydrostatic pressures and re- 
traction speeds, and conducted multiple times to ensure 
repeatability. 

Given the simplicity of the experiment, the resulting 
drop size, r^, shows a surprisingly complex dependence 
on the syringe speed (Fig. [I}. The data can be divided 
into three regions, categorized according the to motion of 
the contact line. At low retraction speeds (Fig. [T]-I), the 
flow from the syringe into the liquid column is relatively 
high, and a bulging liquid bridge forms. The contact an- 
gle on the surface exceeds its equilibrium value, and the 
contact line expands outwards. In this regime, arbitrarily 
large drops can be formed, with the drop radius scaling 
with As one increases the syringe speed (Fig. QJ 

II), the bridge elongation balances the incoming flow and 
the contact line becomes stationary. In this regime, a 
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FIG. 1: (color online) Typical sequence of images of drop 
dispensing at syringe speeds of (I) 15 fim/s, (II) 35 fim/s and 
(III) 400 /j,m/s. The syringe radius is 205 /im. (b) Dispensed 
drop radius, r& vs. syringe speed, U, illustrating the three 
regions: expanding contact line (I), pinned contact line (II) 
and retreating contact line (III). 



quasi-cylindrical liquid bridge forms, finally pinching off 
as it becomes unstable. The resultant drop size still scales 
with U~ x / 2 . However, there is a jump in the drop size 
between the fixed and expanding contact line regimes due 
to a jump in the pinch-off height, h p , (Fig. [5]). A third 
regime is achieved by increasing the syringe speed further 
(Fig. [IJIII) The liquid bridge initially adopts a catenoidal 
shape, but as its contact angle falls below the equilibrium 
value, the contact line retreats inward until, at some criti- 
cal point, the liquid bridge become unstable and pinches 
off. The rapid pinching motion drags the contact line 
inwards at a very high speed (~ 100 times the syringe 
speed), and due to this rapid retreating motion, a small 
drop is deposited, with a diameter that continues to de- 
crease until a minimum drop size is reached, independent 
of retraction speed, with a diameter approximately one 
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FIG. 2: (color online) Pinch-off height, h p , in regions of ex- 
panding (I) and fixed (II) contact line. Numerical simulations 
(□) and experimental measurements (<) are shown. For this 
case, R = 205Aim, Oh = 0.087, Bo = 0.0079. The insets show 
measured and computed shapes of the liquid bridge prior to 
pinch-off (numerical solution: cyan line). 

tenth that of the dispensing syringe. 

We focus first on region II where the contact line is 
stationary, and the bridge near-cylindrical. In this re- 
gion the syringe speed is much smaller than the capillary 
wave speed, U / yf^JpR = O(10~ 4 ), viscous and inertial 
forces are small, and the liquid bridge can be considered 
as quasi-static. Since the Bond number, Bo = gpR 2 /j, is 
small (O(10 -2 )), the hydrostatic pressure in the column 
is approximately constant, dominated by the Laplace 
pressure set by the needle curvature, j/R. This state is 
amenable to the the classical Rayleigh stability analysis, 
which predicts that the height-to-radius, A = h p /R, must 
be less than 2tt to maintain stability [7j. Furthermore, 
for a marginally-stable liquid bridge (i.e. A/2-7T = 1 + e), 
the dimensionless volume, V = v p /irR 2 h p , is given by 
V = 1 + 2e + 5e 2 /2 + 0(e 3 ) ||]. Rearranging this equa- 
tion gives V — (A/27r) 2 + 0(e 2 ), from which we find that 
n p ~ h^. Assuming that the liquid flow from the syringe 
into the liquid bridge, Uf, is constant (driven by the net 
pressure difference, pgH — "f/R), both the volume and 
the height increase linearly with time: v ~ Uft; h ~ Ut, 
from which it is easily derived that the pinch-off height at 
which the bridge becomes unstable, h p , scales like £/ -1 / 2 , 
that the breaking time, t p , scales like C/ -3 / 2 , and that the 
drop radius, r<i, scales like [/~ 1//2 . The experiments sup- 
port this scaling argument very well (Figs. [T]-bj2l and[3]). 

In region I, where the contact line expands, the bridge 
is no longer cylindrical and the Rayleigh stability cri- 
teria can be modified by a small parameter, 8 = (1 — 
w 2 )/(l + w 2 ), where w — r(h]/r(0) is ratio of the upper 
and lower contact line radii [9(. In agreement with the 
theory, we do see a decrease in the pinch-off height, A, 
just as the contact line begins to expand (w < 1), al- 
though S quickly becomes too large for the perturbation 
analysis to remain valid. To address this, we employ a 
numerical model, previously used in studying jet breakup 



3 



10 



V ♦ 



*4 



o 

"D 




10" 



***** 



10"' 10" 
scaled syringe speed U/Uf 



10 



FIG. 3: (color online) Normalized drop sizes, Td/R, at differ- 
ent scaled syringe speeds U/v,f for three syringe radii: 205 /im 
(<), 255 pm (O) and 320 /im (□). The driving pressure, 
Ap = pgH — (7/-R) is 172 (Pa). Also shown are drop sizes 
obtained using an R = 205pm syringe and two effective driv- 
ing pressures: Ap = 52.0 Pa (0) and 291Pa(l>). The syringe 
speed is scaled by the flow speed, Uf = (R 2 /SpL)(pgH —7 / R). 



and the stretching of a pinned liquid bridge [ill with 
modified boundary conditions to include moving contact 
lines. The non-dimensionalized equations for conserva- 
tion of mass and momentum are given by: 



dtr + ur = —ru'/2, 



(d t u 



Wh 



[(r 2 u' 



Bo, 



(1) 



(2) 



while the evolution of the mean curvature, k, (included 
to accurately predict the breakup beyond the validity of 
slender-jet approximation 111 Il9l|) is described by: 



r(l + r' 2 ) 1 /2 (1 + ^)3/2 



(3) 



Here, u{z, t) and r(z,t) are the axial flow speed and 
column radius, normalized by the capillary wave speed, 
u cp = y/j/pR and the syringe radius, R respectively. 
Time, t, is normalized by the capillary time, yj pR 3 /~/. 
A prime denotes the partial derivative with respect to 
the axial coordinate, z. For our experiments, Ohnesorge 
number Oh = fjt/ ^pRy sa O(l). Note that the Weber 
number, We — \/U /u rM , should be small for the equa- 
tions to remain valid [191 ]. 

The boundary conditions require that at the top of 
the liquid bridge, h — Ut, and that the contact line 
is pinned, r(h) = R. The pressure here is the hydro- 
static head minus the pressure drop due to the flow: 
p(h) = pgH — 8pufL/ R 2 (L is the length of the syringe), 
where the inflow velocity, Uf, is evaluated as Uf = u(h). 
At the substrate, z — 0, a solid-wall boundary condi- 
tion is imposed, u(0) = 0. This last boundary condi- 
tion is quite subtle, since the contact line may move with 
time, and this must be consistent with the solution of 
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FIG. 4: (color online) (a) Contact line speed, u c /U, vs. con- 
tact line location, r(0,t)/R, for various scaled syringe speeds: 
U/u f = 0.66(A), 0.99(0), 1.31(t>), 2.63(*), 3.94(v) and 
5.25(0)- The inset shows the corresponding final drop size, 
(b) Corresponding numerical solution. 



the model equations. For a fixed contact line, r(0) is a 
constant and the apparent contact angle varies between 
(9a, 8 r) - the advancing and receding equilibrium con- 
tact angles. For a moving contact line, the contact angle 
deviates from the equilibrium contact angle in order to 
balance the viscous drag. In general this is a function 
of the contact line speed, u c = dtr. Predicting the con- 
tact line behavior and the dynamic contact angle is an 
area of active research (e.g. [la]). For simplicity, we use 
constant values for (6 a ,6 r ) of 100° and 80° respectively. 
The equations {l}-® were solved numerically [u|. Fig. [2] 
indicates excellent agreement achieved between the nu- 
merical solutions and the experimental measurements of 
the pinch-off height, bridge shape, and dependence on 
the syringe speed. 

The speed at which liquid flows into the bridge, Uf, is a 
critical scaling velocity, affecting the pressure boundary 
condition at the top of the bridge and the rate at which 
the bridge volume grows (which in turn plays a central 
role in the bridge stability and subsequent pinch-off). Wc 
can estimate Uf assuming Poiseuille flow through a sy- 
ringe of length L, radius R, and driven by the net pres- 
sure difference, Ap = pgH — "f/R. This estimate differs 
by only a few percent from the value predicted by the 
1-D equations, The appropriateness of this scaling is con- 
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firmed in Fig. [3] which shows the resultant drop size r^/R, 
versus the scaled retraction speed, U/uf, for a series of 
experiments obtained using three syringe sizes and three 
hydrostatic pressures (yielding values of Uf that ranged 
from 46 — 414/um/s). 

The physics of the drop deposition changes abruptly 
in region III where the contact line begins to retreat. 
High speed imaging was used to measure the contact line 
position r(0, i) from which its speed, u c , was calculated 
(Fig. HJa)). Initially, the contact line is at its maximum 
(r(0)/i? w 1), and has zero speed. As the syringe begins 
to retreat, there is a short period of acceleration, after 
which time the contact line moves inward with approx- 
imately constant speed. However, at a critical radius, 
approximately r(Q)/R ~ 0.5, we see a dramatic acceler- 
ation with u c /U reaching 0(100) immediately prior to 
pinch-off. In the constant speed region, u c /u cp -C 1, in- 
ertial and viscous forces are negligible and the contact 
line speed is thus determined solely from the mass bal- 
ance and pressure equilibrium in the bridge. Applying 
this balance, and the fact that our estimate for Uf does 
not depend on U, it is easy to show geometrically that 
u c /U ~ 1 — ut/U, a trend confirmed in both the experi- 
ments and simulations (Fig. [4]). 

The high-speed contact line phase is driven by the 
pinch-off instability, during which time a strong capillary 
force pulls the contact line inward at increasing speed. 
As u c /u cp approaches unity, the viscous forces become 
significant, and the contact line acceleration decreases 
immediately prior to pinch-off. We also see that the ra- 
dial location at which the pinch-off instability initiates 
moves inward as U increases, and that this affects the fi- 
nal drop size (FigH^-inset). A detailed stability analysis 
of the asymmetric liquid bridge with a moving contact 
line explains that the smaller bridge volume correspond- 
ing to the higher syringe speed postpones the pinch-off 
instability, resulting in a smaller drop [20| . However, as 
the syringe speed increases further, the location at which 
the contact line begins its rapid acceleration moves back 
out, and the drop size increases (Fig. Ufa) inset). A pos- 
sible reason for this may be that the dynamic contact an- 
gle decreases, which destabilizes the liquid bridge earlier 
[20l | . The numerical model yields good comparisons with 
the experimental results (Fig. SJd) , capturing the general 
behavior in the constant speed region (including the in- 
crease in u c /U with U) as well as the onset of rapid accel- 
eration prior to pinch-off. However, the model predicts 
the critical radius to be larger than that seen in the ex- 
periment. We believe that the reason for this discrepancy 
is that the numerical simulation uses a static retreating 
contact angle, 9 — 80°, while the contact angle observed 
in the experiment varies, from a value larger than 85° 
prior to the critical point, to a value as low as 60° during 
the acceleration phase. Furthermore, these angles appear 
to depend on the syringe speed, U, and have a strong ef- 
fect on the resultant dynamics [20j. A second reason for 



the degraded agreement between experiment and simula- 
tion may be due to high radial velocities observed which 
violate the assumptions in the the one-dimensional nu- 
merical approach considered here, although despite this, 
the numerical predictions are suprisingly faithful. 

In summary, the retraction speed of the syringe can 
exert a huge influence on the size of the resulting droplet 
that remains on the substrate, with the transition be- 
tween the three regimes identified being determined by 
the balance between the flow into the liquid bridge and 
the onset of the pinch-off instability. At the highest re- 
traction speeds, the small droplet size appears to be de- 
termined by the contact line speed, raising the possibility 
that even smaller droplets might be achieved on surfaces 
that are smoother and/or exhibit higher contact angles. 
The numerical model provides surprisingly accurate pre- 
dictions of the dynamices, even in the regimes where the 
contact line motion and the presence of viscous and iner- 
tial forces make the one-dimensional assumptions ques- 
tionable. 
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